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Abstract 

An analysis method for predicting the inplane stress 
states in anisotropic finite plates with an elliptical 
notch is presented. This method can be used to analyze 
plates with arbitrary notch orientations with respect to 
the plate material axes. The analysis results have been 
validated using finite element analysis results for 
unnotched composite plates and experimental and finite 
element analysis results for stiffened composite panels 
with a skin that has orthotropic properties. The good 
agreement between these results, until the panel 
exhibits nonlinear response either due to bending or 
initiation of damage, indicates that the present analysis 
method can be used to determine accurately the inplane 
stress states and stress concentrations at and around an 
elliptical notch. 

Introduction 

Advanced composite materials are currently being 
considered for aerospace primary structural applications 
due to their potential for reduced structural weight. 
These structures have to be designed to operate safely in 
the presence of structural discontinuities such as holes 
and cracks. Current design and manufacturing practices 
require access holes for assembly and feed-through holes 
for cable and fuel lines. One of the damage tolerance 
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criteria requires that the structure be qualified with a 
crack extending through two stiffener bays. The 
presence of structural discontinuities such as holes or 
cracks can degrade the strength and stiffness of these 
highly-loaded primary structures. Designing efficient 
structures requires a thorough understanding of the 
response of structures with stress concentrations that 
arise due to the presence of holes or cracks and subjected 
to combined inplane loading conditions. 

Analytical solutions to determine stress 

concentrations around holes in anisotropic infinite 
plates were pioneered by Lekhnitski [1] and Savin [2] 
using a complex variable method developed by 
Mushkelishvili. However, no exact closed form 
solutions are available for determining stress 

concentrations in a finite anisotropic plate. Finite 
element analyses are the only methods available for 
obtaining stress distributions in finite composite plates. 
Since finite element analyses are computationally 
intensive, there is a need for developing other analysis 
methods which will enable the designer to evaluate the 
stress concentration problems more rapidly. Applying 
the concept of the Modified Mapping Collocation 
(MMC) technique introduced by Bowie and Neal [3], 
many researchers determined the stress intensity factors 
for cracks in rectangular orthotropic finite plates [e.g., 
4,5] • Sampath and Hulbert [6] determined stress 
distributions in finite orthotropic laminates with 
multiple holes. They also investigated the stress 
concentrations around an inclined elliptical hole in a 
specially-orthotropic laminated plate. Ogonowski [7] 
employed the complex variable technique in conjunction 
with the boundary collocation method to determine the 
stress concentrations around a circular hole in a finite 
laminated plate. Lin and Ko [8] employed the same 
boundary collocation technique to determine the stress 
concentrations and strength predictions for finite 
composite plates with elliptical holes. 

There is no analysis method available which has 
been thoroughly validated using either finite element 
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analysis or experimental results to determine stress 
distributions around an inclined narrow elliptical notch 
in finite anisotropic plates subjected to combined 
inplane loading conditions. The objective of the present 
paper is to present and discuss a validated analysis 
procedure which can be used to determine stress 
concentrations and stress and strain distributions around 
a narrow elliptical notch in finite anisotropic composite 
laminated plates and three-stiffener panels loaded in axial 
compression. The present analysis method was 
validated by comparing the analytical stress and strain 
results with finite element analysis results for square 
laminated plates with notches oriented in two different 
angles to the loading direction. The present analysis 
results and finite element analysis results are also 
compared with the experimental results for three- 
stiffener panels with long elliptical notches. 

Analysis Approach 

This section outlines the governing equations, 
boundary conditions and solution procedure for a finite 
anisotropic plate that has an inclined elliptical notch and 
is subjected to in-plane loads applied along its outer 
edges. The method of analysis is based on (a) classical 
lamination plate theory, (b) anisotropic theory of 
elasticity, and (c) least-squares boundary collocation 
approximation procedure. The geometry and coordinate 
systems used in this study are shown in Fig. 1. The 
plate has a thickness h, length L and width w. The 
analysis method is derived for the general case of an 
elliptical hole having a semi-major axis ‘a’ and semi- 
minor axis ‘b’. The principal axes of the ellipse ate 
oriented at an angle a to the material axes x and y of 
the plate. 

For plane stress problems, the solution is fully 

defined by the stresses O x ,°y, and in a laminated 
composite plate, the in-plane stress resultants and the 
strains are related by 


in which Q a are the reduced stiffnesses and is 
the thickness of the kth layer. For practical 
convenience, the composite laminate is considered as a 
homogenous anisotropic single layer. On combining 
equilibrium and compatibility equations and introducing 
the Airy stress function U(x,y), which is related to the 
stresses as 
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where S# are laminate compliance coefficients and are 

equal to f^ij] . The general expression for the stress 
function U(x,y) depends on the roots of the 
characteristic equation associated with the biharmonic 
equation and is given as 
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Lekhnitski [1] has shown that Equation (5) has no real 
roots. In general, these roots are not equal and can be 
represented by two pairs of complex conjugate roots and 

are given as M-i , M-l , ^2. ^ ^2 • Since the stress 
function is real, the stress function is expressed as 
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Where 
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where U and ^ 2 (Z 2 ) are analytic functions of 
the complex variables Z \ = x + iy and 
Z 2 = X + (X 2 >’ respectively. Introducing these 

( 1 ) functions, the general expressions for stress and 
displacement components can be written as 

°x = 2Re l[x^\(Z]) + q^<P 2 ( Z 2)] 

CT y = 2Re [<p'i( z l) + (p 2< Z 2>] 

(2) T ^ = ± 2Re[Fl(p' 1 (Zi) + F2cp2( z 2)] ^ 
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where 

P l =SllJi? + Sl2 

P 2 = S UP 7 L + S \2 

=Si2l 1 l + — 

M'l 

92 = ^ 12^2 + — 

^2 

The functions q)| (Z | ) and 92 ^ 2 ) are detennined 
by satisfying the boundary conditions. An appropriate 
choice of these functions for a particular problem 
considerably simplifies computational details and speeds 
up the convergence to an accurate solution. The finite 
plate with a hole is considered as a doubly-connected 
region whose outer boundary is rectangular and inner 
boundary is elliptic or circular. For such a doubly 
connected region, the unknown functions can be 
expressed in terms of the Laurent series expansion as 

<Pk(Zk) = A k lnZ k + n g !x B kN zV (*=1-2) 

( 8 ) 

The In Zt term is present whenever the resultant 
of the applied forces on the internal boundary is non- 
zero (i.e., a loaded hole). The coefficients A ^ are 
determined by imposing the conditions of single 

valuedness to the displacements. The coefficients A k 
are equal to zero for an unloaded hole. 

The Laurent series expansion for the unknown 
functions yields a converged solution only for points in 
the complex plane that are a distance greater than or 
equal to a unit radius from the origin. The primary 
complex plane must be conformally mapped onto the 
exterior of a unit circle on a parametric plane to obtain 
converged results. A wide variety of problems can be 
solved once the conformal mapping solution is 
determined for the internal boundary. 

The procedure for determining the mapping 
function is given in Refs. 1 and 2. For a plate with an 
elliptical hole with semi-major and semi-minor axes a 
and b respectively, coinciding with x and y directions, 
the mapping functions are given by 

= z k + \/ z I ± ( a2 + vfo ( k = i , 2 ) 

(a±i\i k b) 



Since the principal axes of the ellipse are oriented at an 
angle a with respect to the principal material directions, 
the parametric equation of the ellipse are 


X = ® cos 0 y' = &sin0 
x = x ' cos a±y' sin a y = x ' sin a + y' cos a 


( 10 a) 

( 10 b) 


where 0 is measured in the x and y coordinate 
system. Substitution of equation (10a) into equation 
( 10 b) results in 

x = a cos a cos 0 ± b sin a sin 0 
y = a sin a cos 0 + b cos a sin 0 .... 


Using the x and y values from equation (11), the 
mapping function given in equation (9) becomes 
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71 1 = (cos 2a + pi sin 2a + Pk sin 2 a) 
n 2 = (sin 2a + Pi cos 2 a ± \i k sin 2 a) 
n 3 = (cos a + sin 2 a) 

^4 = 0 l k cos a ± sin a) 


Equation (12) can be used as the mapping function for 
inclined cracks, where the crack is modeled as an 
elongated ellipse. 

D 

The °kN coefficients of the Laurent series 
expansion are detennined by a least-squares boundary 
collocation technique [Ref. 6 - 8 ], This determination is 
accomplished by limiting the infinite series expansion 
of the stress function to a finite number of positive and 
negative terms. A set of uniformly spaced points is 
selected for the plate and hole boundaries, and the 
boundary conditions are satisfied at these selected 
points. The number of equations for the boundary 
conditions is greater than the number of terms retained 
for the series. A value for the ratio of the number of 
equations to the number of terms of 1.5 to 3.0 is 
suggested in Ref. 6 . The over-determined system of 
equations is represented by the generalized matrix form 


where A denotes the matrix of coefficients from the 
truncated series and B is the vector of boundary 
conditions. The over-determined system of simultaneous 
linear equations is then solved for the unknown 
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coefficients by minimizing the average squared error of 
the boundary conditions. This solution can be 
expressed in a matrix form as 


Results and Discussion 


The results of the present study are presented in two 
A AX = A B (- 14 ) parts. In the first part, the present analysis method is 

where the superscript T denotes the transpose of a validated by comparing the inplane stress results from 

matrix. Once the coefficients of the series are literature and finite element analysis results for 

calculated, stresses can be determined at any point in the finite plates with arbitrarily oriented narrow elliptical 

plate. The truncated Laurent series for the present study notches. In the second part, the present analysis results 

has 23 negative terms and 17 positive terms. Eighty 316 compared with finite element analysis and 

boundary points were selected for the elliptical hole experimental stress results for three-stiffener panels with 

boundary and 20 boundary points were selected for each l° n S anc l narrow notches oriented along two different 

side of the rectangular plate. The choices for the number directions to the loading axis, 
of terms in the truncated series and number of boundary 

collocation points led to convergence for all the results Validation of the present analysis method using finite 
presented herein. element analysis results 


lecimen Description and Test Conditions 


Experimental results from two multi-stiffener 
panels have been used to compare with the present 
analysis results. These specimens were made from the 
Hercules Inc., AS4/3502 graphite-epoxy material 
system. The skin and stiffeners of this panel are made 
from stacks of this material, which were stitched 
together and resin film infiltrated to form a three- 
stiffener panel shown schematically in Fig. 2. Nominal 
values for the mechanical properties for each stack of 
this material are provided in Table 1. Each stack of 
material is 0.058 inches thick and is composed of 7 
plies with a stacking sequence of [45,-45,0,90] s . The 
thicknesses of the plies are 0.00633, 0.01285, and 
0.007018 inches for ±45, 0, and 90-degree orientations, 
respectively. The notch is typically 0.19-in. wide and 
7-in. long and is machined across the panel width 
through the center stiffener either at 90° or at 30° to the 
axial loading direction. Both composite panels are 
stringer-stiffened panels that were machined from the 
upper and lower cover of a 12 -foot-long, stitched 
graphite-epoxy wing box manufactured by the Boeing 
Company (formerly McDonnell Douglas) as part of the 
NASA Advanced Composites Technology (ACT) 
program. A complete description of this wing box is 
given in Ref. 9. The panels were loaded in uniaxial 
compression. The loaded ends of the panels were 
encased in 1.5 inches of epoxy potting compound to 
prevent an end-brooming failure, and the unloaded edges 
of the panel were supported with knife-edge supports to 
provide a simple support boundary condition along 
those edges. A photograph of the test set-up for these 
panels is given in Fig. 3. The electrical resistance 
strain gages used for monitoring the inplane strains are 
also identified in the figure. 


To validate the analysis procedure, a computer code 
was written in FORTRAN 77 and implemented. The 
case of an inclined narrow elliptical hole in a specially 
orthotropic square plate as given in Ref. 6 was solved 
and the results are compared with results from the 
present analysis in Fig. 4. The stress magnitudes 
around the crack radius are normalized with respect to 
the far field stress and these values are plotted as a 
function of 9 in degrees (Fig. 1) to represent the stress 
concentrations. The results obtained using the present 
analysis method compare very well with results in Ref. 
6 . 

A composite plate which is representative of a skin 
on a wing cover panel is considered next. This plate is 
12-in. square and is clamped along all of the edges. This 
plate consists of eight stacks of AS4/3501-6 graphite- 
epoxy material. Nominal values for the mechanical 
properties of this material are provided in Table 1 . Each 
stack of material is 0.058 inches thick and is composed 
of 7 plies with a stacking sequence of [45,-45,0,90^. 
The thicknesses of the plies are 0.00633, 0.01285, and 
0.007018 inches for ±45, 0, and 90 degree orientations, 
respectively. The typical notch is 0.19-in. wide and is 
7-in. long. 

The finite element analysis of the plate is conducted 
using the Structural Analysis of General Shells 
(STAGS) code [10]. The finite element analysis results 
for a case where the notch is inclined at 90° (ct=90°) to 
the loading direction is presented in Fig. 5. The axial 
stress contours shown in Fig. 5(a) indicate that the 
stress concentration is at 6=0°. The axial stress 
magnitudes from the finite element analysis results 
presented in Fig. 5(b) compare very well with the 
present analysis results presented in Fig. 5(c). The 
next case considered has the notch inclined at an angle 
of a=60°. The axial stress results from the finite 
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element analysis presented in Figure 6(a) suggest that 
the stress concentration occurs at 0>O°. The locations 
and magnitudes for the stress concentration predicted 
from the two analyses are shown in Figs. 6(b) and 6(c) 
which indicate good correlation. 

The strain distributions along a line passing 
through the location of axial stress concentration is 
presented in Fig. 7 for notches inclined at a=90° and 
a=60°. These results are plotted as a function of the 
distance from the crack tip, R (see Fig. 1). The axial, 
transverse, and shear strain results from the finite 
element analysis and the present analysis compare very 
well for the range of R values shown in the figure. 
Compared to the case with a=90°, the axial strain 
magnitude decreases and the shear strain magnitude 
increases as the magnitude of the crack orientation angle 
is changed to 60° . This change in stress state in the 
vicinity of the crack tip suggests that the location of 
damage initiation and the damage modes at the crack tip 
could be influenced by changes in the crack orientation 
angle. These results indicate that the present analysis 
method is adequate for predicting the location and 
magnitude of the stresses and strains in composite 
plates with arbitrarily oriented narrow elliptical notches. 

Three-stiffener compression panel results 

Two three-stiffener panels have been tested with 
notch inclination angles of 90° and 60°. A typical panel 
configuration is shown in Fig. 2. A finite element 
model of the test panel is shown in Fig. 8(a) together 
with the finite element discretization in the vicinity of 
the notch tip in Fig. 8(b). To determine the stress state 
in the notch-tip region with the present analysis 
method, the far field axial compressive stress for the 
three-stiffener panel needs to be determined first. When 
the panel has no out-of-plane deformations, the total 
panel load is distributed between the skin and the 
stiffeners in a direct proportion to their inplane 
stiffnesses. The presence of the notch does not 
influence the far-field load distribution in a long panel. 
The far-field axial strains that are calculated using the 
far-field axial stresses are compared with the 
experimentally measured axial strain results in Fig. 9. 
These results suggest that the far-field strain predictions 
compare very well with the experimental results up to 
the point where either the out-of-plane deformations or 
damage initiation at the notch tip cause load 
redistribution. The far-field stresses estimated 
analytically have been used to assess the influence of 
notch angle on the stress distributions around the notch 
tip for the two panels. 


The variation of axial strains at the notch tip for 
the three-stringer panel with a 90° notch is presented in 
Fig. 10 as a function of the far-field stress. The 
experimental results for the axial strain at a point that is 
0.1 -in. away from the notch tip are presented and 
compared with the finite element results and the present 
analysis results. The experimental results suggest that 
a non-linearity exists in the response for far-field 
compressive stress values greater than approximately 
10,000 psi, and this nonlinear response is not captured 
by the finite element analysis results. Since 
geometrically nonlinear finite element analysis was 
performed to predict the panel stress states, this 
departure in the experimental results from the finite 
element analysis results is attributed to the initiation of 
failure at the notch. Damage initiation is not included 
in the finite element analysis. The present analysis 
results differ from the experimental results by a 
maximum of 18 percent until damage initiates in the 
panel. The experimental lateral strain results for this 
test case have not been measured and, hence, are not 
presented here; the shear strain value at the notch tip for 
this panel configuration is zero. 

The experimental and analytical strain results for 
the panel with a 60° notch angle are presented in Fig. 
11. These results are compared at a point that is 0.25- 
in. away from the notch tip. All inplane strain 
components have non-zero values for the case with a 
60° notch orientation. The presence of a combined 
stress state due to the application of a uniaxial 
compression load suggests that the principal stress 
orientations are inclined to the specimen coordinate 
axes. Damage initiation for this test case has been 
observed to occur along a plane that is inclined to the 
mid-section of the panel (Ref. 11). The axial, lateral, 
and shear strain results presented in Fig. 11 (a-c) 
indicate that the panel response becomes nonlinear at a 
far-field stress value of approximately 7,000 psi. The 
finite element analysis results compare very well with 
the experimental results until the far-field stress value 
reaches this level. The axial, lateral, and shear strain 
results from the present analysis are within 10 percent 
of either the experimental or the finite element analysis 
results. The finite element results or the experimental 
results are based on strain values that are smeared over 
an area that is a distance away from the semi-circular tip 
of a slit. The present analysis results are obtained at a 
point in a high strain gradient region along the major 
axis of a narrow elliptical notch. In view of these 
differences, comparisons between the present analysis 
results and the test results are very good. The present 
analysis method is a good preliminary analysis tool to 
predict the stress concentrations around an arbitrarily 
oriented notch. The strain information obtained in this 
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manner can be used to predict damage initiation at a 
notch tip. 

Concluding Remarks 

An analysis approach to determine the inplane 
stress states at the tip of a narrow elliptical notch in an 
anisotropic plate was presented. This analysis approach 
can be used to analyze plates with arbitrarily oriented 
notches and are subjected to combined inplane loading 
conditions. The results from this present analysis 
approach are compared with results from the open 
literature, finite element analysis results presented in 
this paper for laminated composite plates, and finite 
element analysis and experimental results for three- 
stiffener panels loaded in axial compression. The 
present analysis results, finite element analysis results, 
and experimental results are in very good agreement for 
the panels studied until the panels begin to exhibit 
nonlinear response. Geometric nonlinearities and the 
nonlinearities due to damage initiation and growth are 
not included in the present analysis. The correlation 
between the results suggest that the present analysis 
method can be used to determine stress states in 
laminated composite plates with narrow elliptic notches 
and to assess damage initiation for different inplane 
loading conditions. This method also provides a 
verified analysis tool for designing composite structures 
with narrow elliptic and circular cutouts. 
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X 


^Compression load 

A 



R 


a: Angle on inclination for the elliptical notch. 

0: Angular position from the elliptical notch major 
axis. 

R: Distance from the notch tip in inches 


Figure 1 . Plate and notch geometry, dimensions, and 
coordinate system. 



Figure 2. Schematic of the stitched, graphite-epoxy 
panel with a centered, machined notch through the 
center stringer (dimensions are in inches). 
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Figure 3. Photograph of the panel test setup. 
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Figure 4. Comparison of analytical stress results 
around an elliptical notch inclined at 45° to the loading 
direction in a specially orthotropic square plate. 
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Region of high stress 



a. Finite element mesh and location of stress 
concentration 



Region of high stress 

b. Enlarged view of stress contours from finite element 
analysis results 



Region of high stress 
c. Stress contours from the present analysis results 



Region of high stress 

a. Finite element mesh and location of stress 
concentration 

Figure 6. Comparison of axial stress results from finite 
element and present analyses in the vicinity of an 
elliptical notch oriented at a = 60° in an axially 
compressed composite plate. 



Region of high stress 

b. Enlarged view of stress contours from finite 
element analysis results 


Figure 6. Continued. 


Figure 5. Comparison of axial stress results from the 
finite element and present analyses in the vicinity of an 
elliptical notch oriented at a = 90° in an axially 
compressed composite plate. 
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Region of high stress 
c. Stress contours from present analysis results 


Figure 6. Concluded. 
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a. Elliptical notch oriented at a =90° 
Figure 7. Continued. 


Figure 7. Comparison of inplane strain results from 
the finite element analysis and present analysis. 




b. Mesh detail at the notch tip 


Figure 8. Finite element model of the three-stiffener 
panel with a 60° notch. 
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Figure 9. Comparison of experimental and analytical 
far-field axial strain results for the skin. 



Axial compressive strain, p in ./in. 

Figure 10. Comparison of experimental and analytical 
strain results at a point 0.1 -in. away from the notch tip 
for a three-stiffener panel with a 90° notch. 
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a. Axial strain results 


Figure 11. Comparison of experimental and analytical 
axial strain results at a point 0.25 -in. away from the 
notch tip for a three-stiffener panel with a 60° notch. 



Figure 1 1 . Continued. 
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Shear strain, pi in ./in. 


c. Shear strain results 
Figure 1 1 . Concluded. 


11 

American Institute of Aeronautics and Astronautics 



